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Many properties of black holes can be studied using acoustic analogues in the laboratory through 
the propagation of sound waves. We investigate in detail sound wave propagation in a rotating 
acoustic (2 + l)-dimensional black hole, which corresponds to the "draining bathtub" fluid flow. 
We compute the quasinormal mode frequencies of this system and discuss late-time power-law tails. 
Due to the presence of an ergoregion, waves in a rotating acoustic black hole can be superradiantly 
amplified. We compute superradiant reflection coefficients and instability timescales for the acoustic 
black hole bomb, the equivalent of the Press- Teukolsky black hole bomb. Finally we discuss quasi- 
normal modes and late-time tails in a non-rotating canonical acoustic black hole, corresponding to 
an incompressible, spherically symmetric (3 + l)-dimensional fluid flow. 
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I. INTRODUCTION 



Black holes are one of the most fascinating predictions of general relativity. Being a solution of the classical Einstein 
field equations in vacuum, they are the simplest objects which can be built out of spacetime itself. Hawking Q] showed 
that when quantum effects are taken into account, black holes are not really black: they slowly evaporate by emitting 
an almost thermal radiation. Hawking's prediction has been tested time and again in very different ways. It is now 
clear that the appearance of Hawking radiation does not depend on the dynamics of the Einstein equations, but only 
on their kinematical structure, and more specifically on the existence of an apparent horizon 0, The discovery 
of Hawking radiation uncovered a number of fundamental questions: among them the information puzzle, the issue 
of the black hole final state, and so on. For these reasons an experimental verification of the Hawking effect would 
be of the utmost importance. Unfortunately astrophysical black holes, having a temperature much smaller than the 
temperature of the cosmic microwave background, accrete matter more efficiently than they evaporate. An alternative 
could be provided by the existence of large extra dimensions: if gravity is effective at the TeV scale, black holes could 
be produced in particle accelerators However, even if these mini-black holes are actually produced in colliders, it 
is unlikely that they will yield any conclusive evidence of the existence of Hawking radiation. 

Prospects for detecting Hawking radiation changed when Unruh [B| realized that the basic ingredients of Hawking 
radiation can experimentally be reproduced in the laboratory. Since Hawking radiation crucially depends on the 
existence of an apparent horizon, the experimental setup should display the essential features that define apparent 
horizons in general relativity. Unruh considered precisely such a system: a fluid moving with a space-dependent 
velocity, for example water flowing through a nozzle. Where the fluid velocity exceeds the sound velocity we get the 
equivalent of an apparent horizon for sound waves. This is the acoustic analogue of a black hole, or a dumb hole. 
Following on Unruh's dumb hole proposal many different kinds of analogue black holes have been devised, based 
on condensed matter physics, slow light etcetera 0,0- At present the Hawking temperatures associated to these 
analogues are too low to be detectable, but the situation is likely to change in the near future 0, |j| . 

In order to detect Hawking radiation, a full understanding of the classical physics of analogue black holes is 
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necessary. Not only must one control what happens in the experimental situation, but the understanding of classical 
phenomena may bring clues on how to favor the probabilities to detect Hawking radiation. It is also worth stressing 
that some purely classical phenomena shed light on quantum aspects of (analogue and general-relativistic) black hole 
physics. For example, positive and negative norm mixing at the horizon leads to non-trivial Bogoliubov coefficients 
in the calculations of Hawking radiation 10] ; superradiant instabilities of the Kerr metric are related to the quantum 
process of Schwinger pair production Hp : and more speculatively (classical) highly damped black hole oscillations 
could be related to area quantization [12] ] . 

In this paper we carry out a comprehensive study of wave propagation in analogue black holes. We consider in detail 
two acoustic black hole metrics: the "draining bathtub" model for a (2 + l)-dimensional fluid flow and the "canonical" 
metric for a non-rotating, spherically symmetric (3+ l)-dimensional acoustic black hole Q ■ For each acoustic metric we 
compute the characteristic oscillation frequencies using a WKB approach. These complex characteristic frequencies, 
called quasinormal modes (QNMs), play a very important role in the classical physics of black holes. They govern the 
late time behavior of waves propagating outside the black hole: any perturbation of the black hole, after an initial 
transient, will damp exponentially in the so-called "ringdown phase". The frequencies and damping times of this 
ringdown signal depend only on the black hole parameters, such as the mass, charge and angular momentum 13], 
and can therefore be used to estimate black hole parameters from observational data 14]. Therefore, in cases where 
a "formal" definition of (say) the mass and angular momentum of acoustic black holes is missing, quasinormal (QN) 
frequencies may be used to define these quantities operationally. According to recent speculations, highly damped 
QNMs may yield some information on quantum properties of a black hole. In particular, it has been conjectured that 
highly damped QN frequencies could be linked to black hole area quantization [T^. Motivated by these conjectures 
we also study highly damped QNMs of acoustic black holes. We find that for a (2 + l)-dimensional acoustic black hole 
there are no asymptotic QN frequencies, whereas for the (3 + l)-dimensional canonical black hole QN frequencies are 
asymptotically given by Attui = log [(3 — >/E)/2] — i(2n + l)n. This result does not seem to support any of the recent 
conjectures, but perhaps this is no surprise, since Hod's argument relies heavily on black hole thermodynamics. Even 
the very formulation of the laws of black hole thermodynamics for analogue black holes is a non-trivial matter f| . 

After the exponential decay characteristic of the ringdown phase, black hole perturbations decay with a power-law 
tail |l5j ] due to backscattering off the background curvature. Here we compute the late-time tails of wave propagation 
in the draining bathtub and canonical acoustic black hole metrics. We show that for the (2 + l)-dimensional fluid flow 
the field falloff at very late times is of the form $ ~ r' 2m+1 ', where m is the angular quantum number. This time 
exponent is characteristic of any (2 + l)-dimensional flow, and not just of a black hole. For the (3 + l)-dimensional 
canonical acoustic geometry, the field falloff is of the form "J ~ £-(2i+6)_ 

The rotating draining bathtub metric, possessing an ergoregion, can display the phenomenon of superradiance 
|Ttl \yA Hsj . We compute reflection coefficients for this superradiant scattering by numerical integration of the 
relevant equations. Enclosing the acoustic black hole by a reflecting mirror we can exploit superradiance to destabilize 
the system, making an initial perturbation grow exponentially with time: we have an "acoustic black hole bomb" 
[Til I20I I2H . We compute analytically and numerically the frequencies and growing timescales for this instability. 
An interesting feature of acoustic geometries is that the acoustic black hole spin can be varied independently of the 
black hole mass. Therefore, at variance with the Kerr metric, the spin can be made (at least in principle) very large, 
and rotational superradiance in acoustic black holes can be very efficient. It should not be difficult to set up an 
experimental apparatus to observe an acoustic black hole bomb in the lab. A concrete example of an experimental 
setup in which this idea can be realized with moderate experimental effort has been fully worked out by Schiitzhold 
and Unruh [T(|. and is based on the study of gravity waves in shallow water. 

The paper is organized as follows. In Section [H] we discuss the (2 + l)-dimensional draining bathtub metric. We 
first introduce the formalism describing sound propagation in this acoustic metric, and shortly describe the possible 
experimental setup (gravity waves in a shallow basin). Then we compute QNMs, discuss late-time tails, introduce the 
phenomenon of superradiant amplification and quantify the timescales for the acoustic black hole bomb instability. 
In Section ITU we repeat the analysis for the canonical (3 + l)-dimensional acoustic black hole metric (of course, the 
absence of an ergoregion means that we have no superradiance in this case). The conclusions follow in Section ll VI 

II. DRAINING BATHTUB: A ROTATING ACOUSTIC BLACK HOLE 

In this Section we consider a simple "draining bathtub" model, first introduced in [7j, for a rotating acoustic 
black hole. We write down the acoustic metric and the wave equation describing sound propagation in this model, 
and describe a concrete example for a possible experimental setup. We compute QNMs and discuss the late-time 
tail behavior of the system. Due to the presence of an ergoregion, sound waves in this acoustic black hole can be 
superradiantly amplified. We quantify this amplification and discuss the possibility to build an acoustic black hole 
bomb in the lab. 
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A. Formalism and basic equations 



Consider a fluid having (background) density p. Assume the fluid to be locally irrotational (vorticity free) , barotropic 
and inviscid. From the equation of continuity, the radial component of the fluid velocity satisfies pv r ~ 1/r. Irrotation- 
ality implies that the tangential component of the velocity satisfies v s ~ 1/r. By conservation of angular momentum 
we have pv e ~ 1/r, so that the background density of the fluid p is constant. In turn, this means that the background 
pressure p and the speed of sound c are constants. The acoustic metric describing the propagation of sound waves in 
this "draining bathtub" fluid flow is Q : 



ds 2 = - 



c 2 - 



A 2 + B 2 



2A 

dt 2 H drdt - 2Bd(f>dt + dr 2 



(1) 



Here A and B are arbitrary real positive constants related to the radial and angular components of the background 
fluid velocity: 



_ -Ar + B9 

v = . 

r 

This flow velocity can be obtained as the gradient of a velocity potential, v = V?/>, where 

iP = -Alog(r/a) + Bcj). 



(2) 



(3) 



and a is some (irrelevant) length scale. 

In the non-rotating limit B = the metric Q reduces to a standard Painleve-Gullstrand-Lemaitre type metric 
[22 |. The acoustic event horizon is located at m = A/c, and the ergosphere forms at tes — (A 2 + B 2 ) 1 / 2 /c. Unruh 
[5| first realized that the propagation of a sound wave in a barotropic inviscid fluid with irrotational flow is described 
by the Klein-Gordon equation V M V M, J' = for a massless field ^ in a Lorentzian acoustic geometry, which in our case 
takes the form Q. In our acoustic geometry we can separate variables by the substitution 



*(t,r, 0) = i?(r)e i(m ^ 4) . 
Then we obtain a simple ordinary differential equation for the radial variable: 

R rr + Pi {r)Rr + Ql (r)R = , 

where 



(4) 
(5) 



Pi(r) 
Qi(r) 



A 2 + r 2 c 2 + 2\A{Bm - r 2 uj) 
r{r 2 c 2 - A 2 ) ' 
2iABm — B 2 m 2 + c 2 m 2 r 2 + 2Bmuir 2 



(6) 



Notice that if we take the incompressible fluid limit c — > oo we get well known equations from fluid dynamics |23t \24 
We now introduce a tortoise coordinate defined by the condition 



dr* 
dr 



A, 



where A = (1 — A 2 /c 2 r 2 ) 1 . Explicitly, 



A 

2c l ° g 



A 



cr - 



.4 



Setting R = ZH we get: 



ZA 2 H, 



[A(2Z r + PiZ) + A'Z] H, r , + (Z, rr + P\Z,r + Q\Z) H = . 



To obtain a Schrodinger-like equation we impose the coefficient of H tT , to be zero: 



Z., 



(c 2 r 2 - A 2 ) + 2\A(Bm - r 2 w) 
r(c 2 r 2 — A 2 ) 



Z = 0. 



(7) 



(8) 



(9) 



(10) 
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A solution of this equation can easily be found: 

1/9 r / \Bm 
r 1 exp 



1 lnr 



i(A 2 uj — Bmc 2 ) 22 



A 'J"" ' 2Ac 2 
Replacing this solution in Eq. we finally get the wave equation: 

A 2 



ln(cV - A 



1 ( Bm 



1 f 2 1 

3 r -I 



5A 2 



4r 4 c 2 



H = 0. 



(11) 



(12) 



The derivation given here is very similar to the one in [17| . but we have corrected some typos in that paper. In 
particular notice that the potential is the same as equation (10) in [l^, except for a factor 1/4 in the last term. 
Some physical properties of our "draining bathtub" metric are more apparent if we cast the metric in a Kerr-like form 
performing the following coordinate transformation (where again we correct some typos in 0|) : 



dt-^dt- 



Then the effective metric takes the form 



Ai 



-A 2 



dr . 



BA 



r(r 2 c 2 — A 2 ) 



dr . 



ds 2 = - 1 - 



A 2 + £r 



PdP 



1 - 



A 2 



dr 2 - 2Bd<j)dt 



(13) 



(14) 



Notice an important difference between this acoustic metric and the Kerr metric: in the (t, t) component of the 
metric (|14|) the parameters A and B appear as a sum of squares. This means that, at least in principle, there is no 
upper bound for the rotational parameter B in the acoustic black hole metric, contrary to what happens in the Kerr 
geometry. Separating variables by the substitution 



#(t, r, 4>) = % A r i?(r)e i(m *-"* ) 



(15) 



one can show that the radial function H(r) is again a solution of equation l|12|) . 

The wave equation (|12|) can be recast in a more convenient form by the following rescaling: f = rA/c, Cj = ujA/c 2 , 
B = B/A. We get 



QH = 0, 



(16) 



where the generalized potential 



Q={ \ Qj ~^r 



V 



V 



m 

4 



5 

4fl 



(17) 



The rescaling effectively sets A = c = 1 in the original wave equation, and picks units such that the acoustic horizon 
fn = 1- From now on we shall omit hats in all quantities (unless otherwise stated). The rescaled wave equation Ijl6|) 
will be the starting point of our analysis of QNMs, late-time tails and superradiant phenomena. Before giving details 
of this analysis we describe a possible experimental setup in which all of the classical physics we are going to discuss 
could, at least in principle, be reproduced. 



B. A possible experimental setup 

The acoustic metric Q can be realized in a simple experimental setup, that was described in detail by Schiitzhold 
and Unruh [l6l | . In this setup it no longer describes a sonic analogue but rather a shallow basin gravity wave analogue 
of a black hole. The idea is to use gravity waves in a viscosity free, incompressible liquid with irrotational flow: under 
appropriate circumstances, one can envisage the use of common fluids like water or mercury. Schiitzhold and Unruh 
assumed a shallow water, long wavelength approximation: the gravity wave amplitude Sh, their wavelength A and the 
depth of the basin hs are such that Sh -C Kb <C A. Relaxing the assumptions that the bottom of the tank and the 
background flow surfaces are flat and parallel, they showed that the most general rotationally symmetric and locally 
irrotational background flow profile can be described precisely by the draining bathtub metric Q when the radial 
slope of the bottom of the tank - that in cylindrical coordinates (z, r, 4>) will be described by some function f(r) - 
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is small: f'(r) < 1. In this gravity wave black hole analogue the constants A and B are proportional to the radial 
and tangential components of the background flow velocity: 

«* = 4, v T = - h 7== , 2 , (18) 
r/i^/ f + f'{r) 2 

where hoo is the height of the tank far from the black hole, and the slope of the tank satisfies the relation 

/M-JM. (1 „ 

In the previous equations g is the gravitational acceleration, related to the constant c of the acoustic black hole metric 
JU by the relation 



c=VgkZ. (20) 

One of the main advantages of this acoustic black hole is apparent from this equation: the speed of the gravity waves 
can simply be tuned to one's needs by adjusting the height of the basin hoo. Another advantage, that we will not 
exploit here, is that the inclusion of surface tension and viscosity can be used to manipulate the waves' dispersion 
relation. To be concrete we will sometimes consider the following plausible choice of physical parameters, as suggested 
in |l6j |: gravity waves of amplitude Sh ~ 1 mm and wavelength A ~ 10 cm; a tank of height hoo ~ 1 cm (so that 
c ~ 0.31 m s _1 ); and a typical characteristic size of the acoustic black hole horizon r# ~ 1 m, corresponding to 
A = cr H ~ 0.31 m 2 s _1 . 



C. Quasinormal modes 

Numerical and analytical studies of a fairly general class of initial data show that the evolution of the perturbations 
of a black hole spacetime can roughly be divided in three parts: (i) The first part is the prompt response at very 
early times. In this phase, which is the obvious counterpart of light cone propagation, the form of the signal depends 
strongly on the initial conditions, (ii) At intermediate times the signal is dominated by an exponential decay, whose 
frequencies and damping times are determined by the black hole QNMs. In this "ringdown" phase the signal depends 
entirely on the black hole parameters (typically mass, charge and angular momentum), (iii) Due to backscattering off 
the spacetime curvature, at late times the propagating wave leaves a "tail" behind, usually a power law falloff of the 
field. This power law seems to be highly independent of the initial data, and persists even if there is no black hole 
horizon. In fact it only depends on the asymptotic far region. In the following we study the QNM ringing phase of 
our acoustic black hole metric; then we will consider its late-time behavior. 

The characteristic QNMs of the rotating acoustic black hole can be defined in the usual way, imposing appropriate 
boundary conditions and solving the corresponding eigenvalue problem. Close to the event horizon the solutions of 
equation i(Po|) behave as 

H ~ e ±i( ^ Bm)r * . (21) 

Classically, only ingoing waves - that is, waves falling into the black hole - should be present at the horizon. This 
means (according to our conventions on the time dependence of the perturbations) that we must choose the minus 
sign in the exponential. At spatial infinity the solutions of l|16|) behave as 

H ~ e ±iwr * . (22) 

In this case we require that only outgoing waves (waves leaving the domain under study) should be present, and 
correspondingly choose the plus sign in the exponential. This boundary condition at infinity may be cause for 
objections. Indeed, no actual physical apparatus will be accurately described by these boundary conditions: a real 
acoustic black hole experiment will certainly not extend out to infinity. However, we may imagine using some absorbing 
device to simulate the "purely outgoing" wave conditions at infinity (for another example in which an absorbing device 
modeling spatial infinity could be required, cf. Section XI of (16| - in particular their Fig. 5). In any event, later 
on we shall consider the alternative possibility of "boxed" boundary conditions, describing a closed system. Boxed 
boundary conditions were also considered in |2Cj . 

For assigned values of the rotational parameter B and of the angular index m there is a discrete (and infinite) 
set of QN frequencies, loqn, satisfying the wave equation (|l(j|) with boundary conditions specified by Eqs. I|21|) and 
l|22p. The QN frequencies are in general complex numbers, the imaginary part describing the decay or growth of 
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TABLE I: The fundamental (n — 0) QN frequencies for the non-rotating acoustic black hole, using three WKB computational 
schemes. Uq N is the result for the QN frequency using only the lowest approximation |2Sj|. u>q N is the value obtained using 

3rd order improvements pfj]. and finally Wqjv was computed using 6th order corrections |2S||. Notice how for m > 1 the three 
schemes yield very similar answers. 



m 








1 


0.696-0.353i 


0.321-0.389i 


0.427-0.330i 


2 


1.105-0.349i 


0.940-0.353i 


0.945-0.344i 


3 


1.571-0.35U 


1.465-0.353i 


1.468-0.352i 


4 


2.054-0.352i 


1.975-0.353i 


1.976-0.353i 



the perturbation, because the time dependence is given by e~ lbJt . We expect the black hole to be stable against 
small perturbations, and therefore ujqn is expected to have a negative imaginary part, so that the perturbation 
decays exponentially as time goes by. We could not prove stability of the draining bathtub metric in general, but we 
managed to derive upper bounds on the frequencies of unstable modes (if they exist at all). We give the derivation of 
these upper bounds in Appendix lAl As usual, we will order the QN frequencies u)qn according to the absolute value 
of their imaginary part: the fundamental mode (labeled by an integer n — 0) will have the smallest imaginary part 
(in modulus), and so on. 



1. Slowly decaying modes of non-rotating black holes 



The lowest QNMs control the ringing behavior of any classical perturbation outside the black hole. Higher overtones, 
having a larger imaginary part, are damped more quickly and play a negligible role. In particular, it is known |l3| 
that the fundamental (n = 0) QNM effectively determines the response of the black hole to exterior perturbations. 
For these slowly damped QNMs, a WKB approximation - as developed by Schutz, Will and others [2^ l2tt Wh - is 
accurate enough. As a convergence check, we will sometimes extend the WKB treatment to sixth order 28] . The 
WKB method demands that the generalized potential Q defined by Eq. I|17|) uas a single maximum outside the 
horizon. For non-rotating black holes (B = 0) this is certainly true provided m is not zero, and one can easily extract 
the lowest QN frequencies in this case. We show them in Table H] where we also show the convergence of the WKB 
scheme as the order of approximation increases. For m = the situation is not so simple, and we haven't been able 
to extract the QN frequencies using this method, because the potential possesses two extrema. Due to the symmetry 
of the potential l|17fl . QN frequencies of non-rotating black holes are independent of the sign of to (to any positive 
QN frequency for positive m corresponds a negative QN frequency for negative to). 

From Table HI we see that the imaginary part is nearly constant as a function of to, whereas the real part approx- 
imately scales with to. Indeed, in the limit of large to one can show directly from the WKB formula [25j that ujqn 
behaves as 

to 2n + 1 

and indeed already for to = 4 this formula yields very good agreement with the results shown in Table [I] This very 
same result can also be obtained using a Poschl- Teller fitting potential |2^. Just to give an idea of the orders of 
magnitude involved consider an to = 2 mode: if we take the "typical" gravity wave analogue experimental parameters 
of Section Til Bl we get (for the fundamental QNM with to = 2) a frequency lor = 0.945 x c 2 /A = 0.293 Hz and a 
damping timescale r = = 1/0.344 x A/c 2 = 9.37 s. 



2. Slowly decaying modes of rotating black holes 

For rotating black holes, the generalized potential has more than one extremum. This immediately raises problems 
for the applicability of the WKB technique. Previous work on the Kerr geometry j2|| [2^ showed that the WKB 
method can still be used, yielding good results, as long as the rotation parameter is small. The way to handle the 
several extrema of the potential is the following: start with B = and compute the roots of dQ / dr* = (to find the 
extrema; as we said before, for B = the generalized potential Q has a single maximum outside the horizon as long 
as to ^ 0). Compute loqn in the non-rotating case. Then add some small rotation B to the black hole, and compute 
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FIG. 1: In the left panel we show the real part of the fundamental QN frequency to a (and in the right panel we show the 
imaginary part wr) as a function of the rotation parameter B, for selected values of m. For m > (m < 0) ujr and increase 
(decrease) with rotation. 



the roots of dQ/dr* = 0. Now the roots will depend on u>, but only one root yields the correct non-rotating limit as 
B — > 0. It is this root that corresponds to the QNMs of the rotating black hole. Our results, which can be trusted 
for B < 0.2, are shown in Fig. ^ For m > (m < 0) lor and \ujj\ increase (decrease) with rotation, at least in the 
range in which the WKB method can be applied. The QN frequency change with rotation is not dramatic, but it can 
probably be used to apply a "fingerprint analysis" of the acoustic black hole parameters d la Echeverria |l4| : that is, 
once we measure at least two QNM frequencies we may infer the acoustic black hole parameters A and B. Carrying 
out such experiments in the lab can shed some light on the applicability of similar ideas to test the no-hair theorem 
in the astrophysical context [23|. A more complete analysis of the QNMs is still needed to probe the high rotation 
regime: the results shown in Fig. ^seem to indicate that instability may set in for large Bra. A continued fraction 
analysis could be used to test this hypothesis, but it is beyond the scopes of this paper |31j . 



3. Highly damped modes 



QNMs with a large imaginary part, i.e. with a very large overtone number n, have recently become a subject of 
intense scrutiny for black holes in general relativity and similar theories. The interest in these modes comes from 
Hod's proposal [12j that they could be related to black hole area quantization, and from Dreyer's suggestion that a 
similar argument could be used to fix the Barbero-Immirzi parameter in Loop Quantum Gravity j^. In this context, 
an analytical calculation of highly damped QNMs was first carried out by Motl for the Schwarzschild black hole 
[33|. Subsequently Motl and Neitzke [34[ used a complex-integration technique to compute highly-damped QNMs 
of the Schwarzschild and Reissner-Nordstrom black holes. Their analytical results are in striking agreement with 
numerical data [3^, |3|| and alternative analytical calculations [37j ■ The complex integration method has also been 
generalized with success to other black hole geometries j33, H|j , yielding again predictions in excellent agreement with 
the numerical results |40| . In view of this, we shall now build on the results and techniques of [34| to compute highly 
damped QNMs of a rotating acoustic black hole. 

Let us first consider a non-rotating black hole. In the limit B — > the wave equation Hlfi[) reduces to 



d 2 H 

dri 



where the potential 



V(r) 



+ [lu 2 - V(r)] H = 0, 



m 2 - 1/4 5 
j>2 4/* 



(24) 



(25) 



Following (|34|). we can determine the highly damped QNMs looking at the behaviour of the potential near the 
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singular point r = 0. In our case the leading term as r — > is 



^~-Je. (26) 



and in the same limit we have ~ — r 3 /3. Thus 



V — = i- - (27) 

36r2 Arl 1 V ' 

for j = 2/3 (this is precisely the same power-law behavior found in the case of the Schwarzschild black hole, except 
for the different value of j). The result in [34[ carries over directly: 

e 4 ™ = -(l + 2cos7rj). (28) 

However now (1 + 2cos7rj) = 0. This means that there are no asymptotic QN frequencies for this black hole. A 
similar situation occurs for electromagnetic perturbations of a five-dimensional black hole |40j . If we add rotation to 
the black hole, i.e., if B is non-zero, a similar analysis |4l| implies that the asymptotic QN frequencies behave as 

e 47r( w -mB) = _ (1 _|_ 2 COS 7TJ ) , (29) 

where again j — 2/3. Thus also in this case there will be no asymptotic QN frequencies. 

This result is quite puzzling. It could mean either that the real part grows without bound, or that it just doesn't 
converge to any finite value. In any case, an application of Hod's conjecture to these black holes seems impossible: 
the very definition of the classical laws of black hole thermodynamics is non-trivial in the analogue case. In fact, as 
argued by Visser 0, the laws of black hole thermodynamics arise solely from the Einstein equations for the metric; 
from a different perspective, it has been shown that the Einstein equations themselves can be derived from purely 
thermodynamical arguments |4^. A well-known "weakness" of analogue models is the fact that they can be used 
to reproduce the kinematical aspects of the Einstein equations, but not their dynamical aspects (se e eg . |43| for 
the implications of this distinction on the causal structure of acoustic spacetimes). Hod's arguments [l2j assume a 
thermodynamic relation between black hole surface area and entropy: even if Hod's conjecture is true, the absence of 
any such relation for analogues could explain the missing link between the QNM spectrum and area quantization. 



D. Late-time tails 



After the exponential QNM decay characteristic of the ringdown phase black hole perturbations usually decay with 
a power-law tail ^5|, due to the backscattering of waves off the background curvature. The existence of late-time 
tails in black hole spacetimes is by now well established. A strong body of evidence comes from analytical and 
numerical calculations using linear perturbation theory and even non-linear evolutions, for massless or massive fields 
[T5II44I I45ll46l l47| . In a very complete analysis, Ching, Leung, Suen and Young 0,^1 considered the late-time tails 
appearing when one deals with evolution equations of the form <|16|) . and the potential V is of the form 

\rt \ v ( v + 1 ) , cilogr* +c 2 

V{u) - — + , n^oo. (30) 

* * ^* 

By a careful study of the branch cut contribution to the associated Green's function they concluded that in general the 
late-time behavior is dictated by a power-law or by a power-law times a logarithm. The exponents of the power-law 
depend on the leading term at very large spatial distances. The case of interest for us here is when c\ = 0. Their 
conclusions, which we will therefore restrict to the c\ — case, are (see Table 1 in or [IrijV. 

(i) if v is an integer the term v[y + l)/rjj! does not contribute to the late-time tail. Since this term represents just 
the pure centrifugal barrier, characteristic of flat space, one can expect that indeed it does not contribute, at least in 
four-dimensional spacetimes. Therefore, for integer v, it is the C2/r" term that contributes to the late-time tail. In 
this case the authors of [iH lifj find that the tail is given by a power-law, 

* - , /x > 2v + a , a odd integer < 2v + 3. (31) 

where the exponent ji = 2v + 2a — 2. For all other real a, the tail is 



* ~ i-(2f+«) j all other real a 



(32) 
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(ii) if v is not an integer, then the main contribution to the late-time tail comes from the v(y + l)/rj term. In this 
case the tail is 

if ~ ^-(2"+2) j non integer v. (33) 

As an example of non-integer v, Cardoso et al. |47| have shown that the late-time tails of wave propagation appearing 
in odd dimensional spacetimes do not depend on the presence of the black hole at all. For odd dimensions the power- 
law is determined not by the presence of the black hole, but by the very fact that the spacetime is odd dimensional. 
In this case the field decays as 

* „ t-(M+U-2) ; (34) 

where I is the angular index determining the angular dependence of the field, and D the number of spacetime 
dimensions. One can show directly from the flat space Green's function that such a power-law is indeed expected in 
flat, odd dimensional spacetimes. 

From the aforementioned arguments we should expect late-time tails in our draining bathtub acoustic black hole 
to be directly related to the dimensionality of the underlying flow, and not to the presence of a black hole, since we 
are dealing with a (2 + l)-dimcnsional flow. Indeed, at leading order the potential behaves as 

m 2 - 1/4 vjy + 1) 

V - = - , (35) 

^** 

where v = m — 1/2. We are thus in case (ii) above, because v is not an integer. So any perturbation in the vicinities 
of this black hole will die out as a late-time tail of the form 

H ~ t' (2m+1) . (36) 

This power-law is just the one given by the general expression 1)34(1 . if one substitutes D = 3. 

Strictly speaking, the asymptotic form 135|) for the potential is valid only for non-rotating acoustic black holes. For 
the general rotating case the leading term in the potential is w-dependent: 

V ~ m2 I X / 4 = ^ + l) + ^5m 

Since late-times are associated with small frequencies, and the term v(y + 1) gives a contribution to the late-time tail, 
then it follows that the late-time tails of rotating acoustic black holes are the same as those for non-rotating black 
holes, expression l(3f)jl . 



E. Superradiance 

Supcrradiance is a general phenomenon in physics. Incrtial motion superradiance has long been known 48], and 
refers to the possibility that a (possibly electrically neutral) object endowed with internal structure, moving uniformly 
through a medium, may emit photons even when it starts off in its ground state. Some examples of inertial motion 
superradiance include the Cherenkov effect, the Landau criterion for disappearance of superfluidity, and Mach shocks 
for solid objects travelling through a fluid (cf. |49lfor a discussion). Non-inertial rotational motion also produces 
superradiance. This was discovered by Zel'dovich |5Cj, who pointed out that a cylinder made of absorbing material 
and rotating around its axis with frequency f2 can amplify modes of scalar or electromagnetic radiation of frequency 
u), provided the condition 

lu < mfi (38) 

(where m is the azimuthal quantum number with respect to the axis of rotation) is satisfied. Zel'dovich realized 
that, accounting for quantum effects, the rotating object should emit spontaneously in this superradiant regime. He 
then suggested that a Kerr black hole whose angular velocity at the horizon is O will show both amplification and 
spontaneous emission when the condition l|38|) for superradiance is satisfied. This suggestion was put on firmer ground 
by a substantial body of work |5l| . In particular, it became clear that (even at the purely classical level) superradiance 
is required to satisfy Hawking's area theorem |19l l52| . 

Superradiance is essentially related to the presence of an ergosphere, allowing the extraction of rotational energy 
from a black hole through a wave equivalent of the Penrose process 0. Under certain conditions, superradiance can 
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be used to induce instabilities in Kerr black holes 11]. Indeed, all spacetimes admitting an ergosphere and no horizon 
are unstable due to rotational superradiance. This was shown rigorously in |54|. but the growth rate of the instability 
is too slow to observe it in an astrophysical context 153. Kerr black holes are stable, but if enclosed by a reflecting 
mirror they can become unstable due to superradiance |2oL l5c| ; we will discuss this "black hole bomb" instability as 
applied to the analogue acoustic black hole. 

The possibility to observe rotational superradiance in analogue black holes was considered by Schutzhold and Unruh 
fl(| . and more extensively by Basak and Majumdar 0,^3' wno computed analytically the reflection coefficients in 
the low frequency limit uiA/c 2 <C 1. In particular, the authors of [l6l ] showed that the ergoregion instability in gravity 
wave analogues is related to the existence of an "energy function" [their Eq. (68)] that is not positive definite inside 
the ergosphere. In the context of analogues, inertial superradiance based on superfluid 3 He has been studied by 
Jacobson and Volovik pS7|. 

Here we present a quantitative calculation of the efficiency of superradiant amplification for the draining bathtub 
metric. We work in the frequency representation, and consider an incident plane wave of unit amplitude at infinity, 
frequency to and azimuthal index m. Part of this wave will be reflected back by the medium, the reflection coefficient 
being some (complex) number R um - In terms of the wave equation (|16J) . this determines the following boundary 
condition at infinity: 



At the sonic horizon (r — ► 1, r* 



-oo ) the solution behaves like 



-i(w-mS)r* 



l 



(39) 



(40) 



where T^ m is the transmission coefficient. An easy way to prove the existence of superresonance is to compute the 
Wronskian of a solution of equation (|16[) and of its adjoint at the sonic horizon and at infinity. From the constancy 
of the Wronskian as a function of the radial coordinate, using the boundary conditions i|39[l and (|40|l we find the 
following "energy conservation" condition: 



/ mB 

1 - \RumA = ( 1 

UJ 



IT, 



(41) 



Therefore, if u> < mB the reflection coefficient |-R W ro| > 1, i.e. we have superradiance. 

To compute numerically with improved accuracy the reflection coefficient i? wm , we have used a refined condition 
at the horizon: 



-i(u-mB)r t [1 +yi ( r _1) + 



(42) 



where the leading-order coefficient is 



(1 + 2B 2 )m 2 - 2ujBm + 1 
2[i(mB - oj) + 1] 



(43) 



We also keep higher-order terms to extract the reflection coefficient at infinity. Namely, for the outgoing wave we use 
an expansion of the form: 



H ~e u 



1 



zi , z 2 



(44) 



where 



Zl = 1- 



,5 + (4m 2 - 1) 
~8u 



z 2 = IZi 



Bm (4m 2 - 9) 



16w 



and for the ingoing wave we use the complex conjugate of Eq. (|44|l . As a check of our code we have reproduced results 
by Andersson et al. |58j ] for the superradiant amplification of a scalar field in the vicinities of a Kerr black hole. Our 
results are in perfect agreement with Fig. 1 of |5Sj . In particular, we find a maximum amplification coefficient given 
by 1 — |i? wm | 2 ~ 0.2 % for I = m = 2 scalar perturbations of a near-extremal extremal Kerr black hole (in early 
numerical work |l9ll56j] the maximum amplification was found to be 1 — |i? wm | 2 ~ 0.3 %). 

Results of the numerical integrations for the draining bathtub metric are shown in Fig. Panels on the left show 
the reflection coefficient |i? w i| 2 for m = 1, and panels on the right show |i? U 2| 2 for m = 2, for selected values of the 
black hole rotation B. Panels on top show that, as expected, in the superradiant regime < to < mB the reflection 
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FIG. 2: Reflection coefficient \Rujm\ 2 as a function of uj for m = 1 (left panels) and m = 2 (right panels). Each curve 
corresponds to a different value of B, as indicated. The top panels show that the reflection coefficient decays exponentially at 
the critical frequency for superradiance, ujsr = mB. The middle panels show a close-up view in the superradiant regime for 
B < 1: at B = 1 the maximum amplification is 21.2 % (m = 1) and 4.7 % (m = 2). The bottom panels show that superradiant 
amplification can become much more efficient for values of the rotation parameter B > 1. 
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coefficient |-Ru, m | 2 > 1. Furthermore, as one increases B the reflection coefficient increases, and for fixed B, the 
reflection coefficient attains a maximum at w ~ mB, after which it decays exponentially as a function of u 

outside the superradiant interval. This is very similar to what happens when one deals with massless fields in the 
vicinities of rotating Kerr black holes |l9j| . In particular, from the close-up view in the middle panels we see that, for 
B = 1, the maximum amplification is 21.2 % (to = 1) and 4.7 % (to = 2). 

As a final remark, and as we have anticipated, an important difference between the acoustic black hole metric and 
the Kerr metric is that in the present case there is no mathematical upper limit on the black hole's rotational velocity 
B. In the bottom panels we show that, considering values of B > 1, we can indeed have larger amplification factors 
for acoustic black holes. 

Summarizing: if we are clever enough to build in the lab an acoustic black hole that spins very rapidly, rotational 
superradiance can be particularly efficient in analogues. This is an important result, considering that the detection of 
rotational superradiance in the lab is by no means an easy task, as originally predicted by Zel'dovich and confirmed 
by recent reconsiderations of the problem [49^ . Of course, in any real-world experiment the maximum rotational 
parameter will be limited. At the mathematical level, the equations describing sound propagation (which are written 
assuming the hydrodynamic approximation) will eventually break down. Physically, if the angular component of the 
velocity v e becomes very large the dispersion relation for the fluid will change, invalidating the assumptions under 
which we have derived our acoustic metric |l6j |. To be more concrete, let us consider again the gravity wave analogue 
described in Section lll Bl Then we can use a very simple argument to limit the acoustic black hole spin. The derivation 
of the rotating acoustic black hole metric is based on the assumption that the bottom of the tank should not be too 
steep, that is, f'(r) <C 1. This condition translates into a condition for B: 

/ 3 \ 1/2 

B « - A 2 J , (45) 

which must be satisfied for all values of r, and in particular at the acoustic black hole horizon r = Th (notice that 
here, and here only, we have switched back to physical units). Setting r = rn in the previous inequality and using 
the physical parameters quoted in Section Til Bl we get the very stringent condition that B < 2.2 m 2 s _1 , or (in the 
dimensionlcss units we use throughout this paper) B <C 7.0 . In other words, we can only get values of the rotation 
parameter larger than B ~ 1 when the slope of the tank is so large that the assumptions underlying the derivation 
of the acoustic metric are not valid any more. Of course, this example does not mean that such a constraint applies 
to every possible experimental realization of the draining bathtub metric. However, it serves as an illustration of the 
kind of experimental difficulties we may expect to encounter in practice. 

The superradiant phenomena we have described are purely classical in nature. However, an interesting suggestion 
to observe quantum effects in acoustic superradiance was put forward in |l8j . To write down our acoustic metric we 
required the flow to be irrotational and nonviscous. As a natural choice, we could use a fluid which is well known 
to possess precisely these properties: superfluid Hell. In this case the presence of vortices with quantized angular 
momenta may lead to a quantized energy flux. The heuristic argument presented in goes as follows. Let us 
imagine that our black hole is a vortex with a sink at the centre. In the quantum theory of Hell the wavefunction 

is of the form ^ = exp i^ </>(rj- Aground , where fj is the position of the j-th particle of Hell. The velocity at any 

point is given by the gradient of the phase at that point, v = V(j), so that (roughly speaking) the velocity potential 
(0 can be identified with the phase of the wavefunction. This phase will be singular at the sink r = 0. Continuity of 
the phase around a circle surrounding the sink requires that the change of the wavefunction satisfies A<j> = 2-kB. For 
the wavefunction to be single valued, B (that is, the black hole's angular velocity at the horizon) must be the integer 
multiple of some minimum value AB, i.e., B = n AB. Then the angular momentum of the acoustic black hole would 
be forced to change in integer multiples of AB. Correspondingly, the spectrum of the reflection coefficients may be 
given by equally-spaced peaks with different strengths. This discrete amplification could enhance chances of observing 
superradiance in acoustic black holes, and rule out (or provide empirical support to) some of the many competing 
heuristic approaches to black hole quantization. 



F. Superradiant instabilities: the acoustic black hole bomb 



The very existence of an ergoregion in the acoustic black hole metric allows immediately for the possibility to make 
the system unstable. Suppose we enclose the system by a reflecting mirror at constant radius tq- Now, throw in a 
wave having frequency ui — ujr + iwj such that lor < mQ: the wave will be amplified at the expense of the black 
hole's rotational energy and travel back to the mirror. There it will be reflected and move again towards the black 
hole, this time with increased amplitude. Through repeated reflections, the waves' amplitude will grow exponentially 
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FIG. 3: Real part (left) and imaginary part (right) of the fundamental BQN frequency and the first two overtones as a function 
of mirror location ro. Both panels refer to m — 1 and B — 1 (but the B-dependence of the real part of the BQN frequency is 
very weak). With excellent accuracy u>i crosses zero (and the instability abruptly shuts off) at the critical radius predicted by 
Eq. <EHJ, as can be verified taking a glance at the first row of Table El 



with time. This "black hole bomb" was first proposed by Press and Teukolsky |56(, and recently it has been studied 
in detail by some of us [2(j in the context of the Kerr metric. The detailed analysis performed in [2(j showed that 
the black hole bomb can be characterized by a set of complex resonant frequencies, the Boxed QuasiNormal Modes 
(BQNMs). The real part of a BQNM is, not surprisingly, proportional to 1/Vo, where r is the mirror radius: this is 
essentially the condition for the existence of standing waves in the region enclosed by the mirror. The imaginary part 
of the BQNMs is proportional to (uir — mfi): the system can only become unstable if it is in the superradiant regime. 
Combined with the standing wave condition this implies that the mirror should be placed at some radius ro > 1 /mil 
in order for the system to be unstable. As rotational energy is extracted from the system the black hole will spin 
down. For any given ro the instability will eventually shut off, as the condition ro > 1/mfi will no longer be satisfied. 

An additional reason to study the system "acoustic black hole+mirror" is related to the problem of boundary 
conditions at infinity for QNMs (cf. Section lil Cjl . Whatever analogue model we consider, no physical apparatus will 
ever extend out to infinity. We are left with two possibilities: we can either use some absorbing device to simulate 
spatial infinity, or we can impose alternative boundary conditions on the system. A natural choice is to have a 
reflecting mirror surrounding the apparatus, so that we are effectively building an acoustic black hole bomb. 

In the following we will assume that we are in the presence of a perfectly reflecting mirror. Correspondingly, we 
shall impose the boundary condition H = at r = ro . At the horizon we demand, as usual, the presence of purely 
ingoing waves (waves headed towards the horizon), i.e. H ~ e -i(w-Bm)r, _ -y^ e are a g am m presence of an eigenvalue 
problem. However, since the boundary conditions have changed we shall not refer to the characteristic frequencies 
as QN frequencies, but rather as Boxed QuasiNormal frequencies (BQN frequencies). The reader is referred to [20|. 
where this terminology was first introduced. Using a direct integration of the wave equation (see Section III El for 
details) it is quite easy to compute the BQN frequencies. We simply integrate Eq. (|16fl outwards from the horizon, 
where we impose the condition (|42|l . to the mirror location ro. The (complex) BQNM frequencies are those frequencies 
for which 



H(w B qnm, r ) = . (46) 

Some results for the fundamental BQNM and the first two overtones are shown in Fig. [3] There we fixed the black 
hole rotation parameter B = 1, and m = 1, as a typical case. 

From the left panel we see that the real part of the BQNM frequency for the n-th overtone (where we use the 
convention that n = corresponds to the fundamental mode) scales as 1/Vq, as anticipated. The proportionality 
constant can be obtained by analytical arguments similar to those in |20| . The result is 

Jm,n / a i-t\ 

ur = , (47) 

where the j m ,nS are zeros of the Bessel function of integer order m (recall that for a Kerr black hole bomb one gets 
= ji+i/2,n/ r o instead, cf. |2£})- The numerical values of the j m ,n' s can be found, e.g., in [H^. For reference, we list 
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TABLE II: Zeros of the Bessel functions j m ,n for the first few values of m. The real part of the BQNM frequency is well 
approximated by Eq. 14711 . 



m 


jm,a 


jm,l 


Jm,2 


1 


3.83171 


7.01559 


10.17347 


2 


5.13562 


8.41724 


11.61984 


3 


6.38016 


9.76102 


13.01520 


4 


7.58834 


11.06471 


14.37254 


5 


8.77148 


12.33860 


15.70017 




the first few values in Table ITU Our numerical calculations are in excellent agreement with the predictions of formula 
(|47|l . This is particularly true for small values of B, but the -B-dependence of the real part of the BQN frequencies is 
very weak anyway, as it is for the Kerr metric pflj . 

From the right panel we see the behavior predicted by the qualitative arguments at the beginning of this Section: 
for any BQNM, the instability is more and more efficient as the mirror radius r$ becomes smaller, until eventually 
the mirror radius becomes small enough that the instability shuts off. 

A few remarks are in order: (i) The real part of the BQN frequency (slowly) increases with overtone number n. 
(ii) Higher overtones become stable at larger distances, but they also attain a smaller maximum growing rate, (iii) 
With excellent accuracy, the instability switches off at the critical radius predicted by the analytical formula 147(1 
supplemented by the superradiance condition (|38fl. that is: 

ro,c - ■ (48) 

Since j m ,n grows linearly with to (in particular, this is asymptotically true for high to's i59j) the critical radius is 
almost to- independent. However, the instability is not so efficient for higher to as it is for small m, at least when 
B < 1. This is apparent in Fig. 0] where we compare growth timescales with to = 1 (left) and m = 2 (right) for 
acoustic black holes rotating at different rates. For example, the maximum growth rate for B = 0.3 is ~ 10~ 4 for 
to = 1, and ~ 10~ 5 for to — 2. The growth rate for m = 2 becomes roughly comparable to the growth rate for m = 1 
when B > 1: look for example at the curves corresponding to B = 5.0. However, other physical considerations may 
forbid the construction of acoustic black holes rotating at such rates (see the arguments at the end of Section Hi Ell . 

To get an idea of the orders of magnitude involved, let us consider again the "typical" parameters for the gravity 
wave analogue we introduced in Section [IIBI Let us pick m = 1 and an acoustic black hole rotation parameter B = 1, 
which - according to the arguments in Section III El - is close to the maximum rotation rate we may hope to achieve 
(recall hats distinguish dimensionless quantities from quantities in physical units). To take advantage of the process 
the mirror should be located close to the maximum of w/, but not quite at the maximum. For example, if we place 
the mirror at rg ~ 10 (cf. the right panel of Fig. [3J we get Cjj ~ 4 • 10~ 3 (in this case the maximum growth rate 
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would be <2jj — 8.5 • 10~ 3 ). This corresponds to a growth time r ~ 800s ~ 13 minutes (at the maximum we would 
get r ~ 6 minutes). Roughly speaking, this means that the amplitude will double every 13 minutes or so, and will be 
amplified by orders of magnitude on timescales of the order of a few hours. To reduce the typical timescale to seconds, 
one just has to ajust the horizon radius and wave velocity. For example, working with rn ~ 0.1 m and c~10ms , 
we would get a typical timescale of about 2 seconds. This looks like a perfectly reasonable timescale to observe an 
acoustic black hole bomb in the lab. 

When the instability sets in the acoustic black hole loses energy and angular momentum (AE ~ BAJ). The 
critical radius r$ tC increases with decreasing B. To estimate the efficiency we can use a simple argument. We can 
take advantage of the mechanism by picking some large B and ro ~ 2ro^ max , where ro tmax is the radius of maximum 
growth timescale at the given B. As the hole emits angular momentum B will decrease. For any given mirror location 
ro, the bomb eventually switches off when condition 148|) is satisfied. From the difference in initial and final angular 
momentum AB we can infer the extracted energy AE. For more details we refer to pcj. 

As a final remark we want to stress that, when we talk about surrounding the acoustic black hole by a reflecting 
mirror, what we have in mind is a generic type of mirror. For instance, in the shallow basin gravity wave analogue, 
the mirror could be a circular rubber band of radius ro reflecting gravity waves. A possible alternative to implement a 
mirror could involve, for example, variations in the dispersion relation for sound waves. In particular, Schiitzhold and 
Unruh |l6j | suggested that a simple change of the height hoc of the basin would imply a change of the gravity wave 
speed, hence of the effective refractive index. The changed refractive index may be used to trap the waves, which can 
then be amplified via superradiant scattering until non-linear effects dominate. 



III. THE CANONICAL NON-ROTATING ACOUSTIC BLACK HOLE 



In the previous Section we considered the metric describing a rotating (2 + l)-dimensional acoustic black hole. 
Rotation implies the presence of an ergosphcrc, allowing us to study in the lab some of the most interesting phenomena 
concerning black hole physics (eg. superradiance and the related rotational instabilities). However, the non-rotating 
limit of the acoustic metric we considered is not the "natural" metric describing a non-rotating acoustic black hole. 
In this Section we are going to consider another class of non-rotating, (3 + l)-dimensional acoustic metrics. These 
metrics are associated to the most general spherically symmetric flow of an incompressible fluid: in this sense they 
represent "canonical" metrics for a non-rotating acoustic black hole. 



A. Formalism and basic equations 



It is simple to show that the acoustic metric corresponding to the most general spherically symmetric flow of an 
incompressible fluid is given by 



ds A 



dt 2 



dr 



2 i 2 



sinr 



(49) 



This metric does not correspond to any of the geometries typically considered in general relativity, but it describes 
(in the sense specified above) a "canonical" acoustic black hole. The propagation of small disturbances (sound waves) 
is again described by the massless Klein-Gordon equation V A1 V M, I , = in this background. We can separate variables 
by the substitution 



$>(uj r) ■ + 
*(t,r,<t>) = ^^-e-™ t Y lm (9), 



(50) 



where Yi m (9) are the usual spherical harmonics. This yields the wave equation 

,2 



— - V $ = 



where 



V 



1 



1(1 + 1) 



4r 4 



(51) 



(52) 



and the tortoise coordinate r* is defined, as usual, by dr/dr* — (1 — r^/r 4 ). Since this is a spherical symmetric 
problem, the azimuthal number m does not play any role, only the angular momentum I is important here. In the 
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TABLE III: The fundamental (n = 0) QN frequencies for the non-rotating canonical acoustic black hole, using three WKB 
computational schemes. WqJ, is the result for the QN frequency using only the lowest approximation |25j|. Uq], is the value 
obtained using 3rd order improvements pfj . and finally ujqI^ was computed using 6th order corrections psj . 



/ 


^QN 









1.13-0.73i 


0.19-1. Hi 


0.06+0.87i 


1 


1.37-0.68i 


0.53-0.7H 


1.09-0.39i 


2 


1.82-0.64i 


1.41-0.6H 


1.41-0. 70i 


3 


2.36-0.62i 


2.10-0.62i 


2.12-0.62i 


4 


2.94-0.62i 


2.75-0.62i 


2.75-0.62i 



next Section we summarize our results for QNMs and wave tails of the canonical acoustic black hole. When presenting 
our numerical results we will choose units such that c = tq = 1. This is of course equivalent to a simple rescaling of 
the radial variable (f = r/ro) and of the frequency (Q = iv/c), but we will omit hats in the following. 



B. Quasinormal modes 

A 6th order WKB analysis reveals an unstable QNM for I = 0. However, it is possible to prove stability of the 
canonical black hole since the potential is positive-definite. We therefore computed the QN frequencies using first the 
lowest approximation |25|. then 3rd order corrections |26| and finally 6th order corrections. The results are presented 
in Table ITTT1 The QN frequencies for I = and I = 1 seem to be the problem here. For any other I the value quickly 
converges as we increase the correction order of the WKB method. Notice for example that the I = mode suffers a 
variation of almost an order of magnitude as we go from the lowest approximation to the 3rd order correction scheme, 
and gets unstable for the scheme using 6th order corrections. The reason for this failure is most likely related to the 
breakdown of the basic WKB assumptions, namely that the ratio of the derivatives of the potential to the potential 
itself is small. Indeed a close inspection shows that as I increases these ratios tend to decrease. This means that the 
WKB method is more reliable for higher I: this was first observed in the early works on the subject |2g,|2g. For I = 
and I = 1 the lowest WKB approximation gives the most reliable results. This is confirmed by recent numerical work 

In the limit of large I one finds 

lV2 . yg(l + 2n) 
W ~ ^37? ~ 1 3 3/4 > ( 53 ) 

a result which agrees very well with our WKB data already for 1 = 3. 

The calculation of highly damped QNMs proceeds along the same lines sketched for the (2 + l)-dimensional acoustic 
black hole. In this case the index j — 3/5, which implies that 

3 — \f% 

Attu = log — i(2n + l)ir . (54) 

Once again asymptotic QN frequencies are not given by the logarithm of an integer, as required in Hod's construction 
[l2|. As we remarked earlier this does not necessarily imply that the conjecture is wrong, since a thermodynamical 
interpretation of the black hole area is possible only when the dynamics of the system are described by the Einstein 
equations. 



C. Late-time tails 



The analysis of late-time tails proceeds as in the case of the (2 + l)-dimensional acoustic black hole. We can easily 
show that asymptotically the potential behaves as 



(55) 
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Following the previous analysis we now find a = 6, and since / is an integer the power-law falloff is of the form 

$ ~ t-( 2/+6 ) . (56) 

Thus any perturbation eventually dies off as £~( 2J + 6 ) 5 much more quickly (for m = I) than in the (2 + l)-dimensional 
acoustic black hole background. 

IV. CONCLUSIONS 

In this paper we have considered two acoustic black hole metrics: the (2 + l)-dimensional "draining bathtub" metric, 
and the (3 + l)-dimensional canonical non-rotating acoustic black hole. We have studied QNMs and late-time tails in 
both metrics, and superradiance in the "draining bathtub" metric. 

More specifically, we numerically computed slowly damped QNMs of the draining bathtub metric using a third- 
order WKB approach for all values of m ^ 0. We analytically evaluated QN frequencies in the large-m limit [Eq. 
(EH] and set upper limits on frequencies of unstable modes f Appendix IX}) . We showed that highly damped modes do 
not tend to any simple limit. At late times, power-law tails decay as t~( 2m+1 \ This behavior is typical of any odd- 
dimensional spacetime: independently of the presence of a black hole, if D is odd the power-law falloff is proportional 
to t-& l+D -V [the (2 + l)-dimensional case, D = 3, is special: the azimuthal and angular numbers are the 
same, I — m, since there is only one angular coordinate]. The draining bathtub metric, possessing an ergoregion, 
can superradiantly amplify waves. We computed reflection coefficients for this superradiant scattering by numerical 
integration. When the (dimensionless) acoustic black hole rotation B = 1, the maximum amplification is 21.2 % for 
m = 1 and 4.7 % for m = 2. Enclosing the acoustic black hole by a reflecting mirror we can destabilize the system, 
making an initial perturbation grow exponentially with time: we have an acoustic black hole bomb jlflj . or "dumb 
hole bomb" . We computed analytically and numerically the frequencies and growing timescales for this instability. 
An interesting feature of acoustic geometries is that the acoustic black hole spin can be varied independently of the 
black hole mass. Therefore, at variance with the Kerr metric, the spin can be made (at least in principle) arbitra rily 
large, and rotational superradiance in acoustic black holes can be very efficient. Gravity waves in shallow water [l(j 
provide a concrete example of an experimental setup for studying classical physics in a draining bathtub metric [and 
of experimental difficulties in increasing arbitrarily the rotation rate: see the discussion leading to Eq. I|45|l ]. In this 
case, and for a typical choice of parameters, growth times for the black hole bomb instability would be of the order of 
minutes: this appears to be well within the range of experimental possibilities. Finally one can speculate, following 
[l8|. that superfluid Hell could be used to observe some quantum effects, such as a discrete superradiant amplification. 

For the canonical acoustic black hole we computed slowly-damped QNMs using first, third and sixth order WKB. 
For I < 2 the WKB method does not seem to converge, and even yields an unstable frequency at 6th order when 
I = 0. This is only due to bad convergence properties of the WKB technique, since the canonical acoustic black 
hole is stable. In the high-damping limit QNMs are given by Attlu = log [(3 — \/5)/2] — i(2n + 1)tt, so they are not 
proportional to the logarithm of an integer, as required by recent conjectures. This is not surprising since there are 
no Einstein equations for acoustic metrics, i.e., the acoustic metric evolution is not governed by Einstein's equations. 
Therefore, acoustic black holes can teach us a lot about quantum gravity but cannot shed any light on its 

possible connection with highly damped quasinormal modes. 

The late time falloff in the canonical acoustic black hole metric is proportional to £-( 2Z + 6 ) ; to be compared with 
the t~( 2l + 3 ) decay of 4-dimensional Schwarzschild black holes and with the £-( 2 '+ 3D -s) decay of even-dimensional 
spherically symmetric black holes with D > 4 |47T |. 
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APPENDIX A: UPPER BOUNDS FOR QUASINORMAL FREQUENCIES OF UNSTABLE MODES 

To derive some bounds on the magnitude of QN frequencies of possible unstable QNMs we shall follow Detweiler 
and Ipser jgjj. Let us begin with the Klein-Gordon equation for the evolution of the field, V /J V 7i \I r = 0. Using the 
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metric (|14|) and performing a mode decomposition \&(t, r, 0) = R{r)e^ m ^ ut \ this can be written as 

r 9 d , „ n 2Bmuj / 1 B 2 \ 9 , A . 

—u> 2 R+- [rfB! ]--—R-( \ m 2 R = 0j Al 

jcr dr c z fr \r c z r 6 j J 

where we have defined / = A -1 = 1 — A 2 /c 2 r 2 . 

Multiply by R* and integrate from the horizon to spatial infinity. The result is 



dr 



-V|i?| 2 - rf\Rf - 2 -^p\Rf - (--^s?) '»'-' I*!'" 
jc c jr \r err j 



0, (A2) 



where we have used an integration by parts and discarded the surface integrals. In fact, for unstable modes (and for 
these only) the boundary conditions guarantee that an unstable mode vanishes exponentially as r — ► r# or r — ► oo. 
The imaginary part of this equation yields 

dr (r 2 u R - Bm) ^ = , (A3) 



frc 2 

where we used u> — lur + iuji. Therefore ujr and Bm must have the same sign for unstable modes. Furthermore, since 
\R\ 2 / frc 2 is always positive, the quantity {r 2 u>R — Bm) must be negative somewhere for the integral to vanish. Since 
r 2 u)R increases with r, if this quantity is somewhere negative, then it certainly is negative at the horizon. In this way 
we get an upper bound for ojr: 

^B<-2-- (A4) 

' H 

To get a similar bound on wj consider now the real part of Eq. I|A2|I : 



dr 

fff 



2Bmw R c 2 fm 2 B 2 m 2 9 Ji?'| 2 l r\R\ 2 

fc 2 



R 



= . (A5) 



The positive part of the integrand is (assuming without loss of generality that Bmcuji > 0) 

2 2Bmoj R c 2 fm 2 ,„ 2 |i?'| 2 2 

"i + + -75- + / c > ^ ■ (A6) 

The negative part is maximized at the horizon: 

9 B 2 m 2 9 B 2 m 2 , , , 

w « + —r- < + —r- ■ ( A? ) 

' ' H 

A necessary condition for the integral to vanish is that the integrand be negative at the horizon, and this implies 

9 , B 2 m 2 . 

, A8 

r H 

Using the bound (|A4|) for lur we finally get an upper bound for u>i: 

9 2B 2 m 2 

lj 2 < . A9 

Notice that the upper bound for ujr ensures we are in the superradiant (or better, superresonant) regime. This seems 
to be a general feature, since it has been shown to hold in several black hole spacetimes. If the instability sets in at 
all, the system behaves as a black hole bomb. 
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